o 

(N 

C 
=3 



cr. 

bJQ, 



> 
in 
in 

in 
d 

On 
O 



X 



IFT-UAM/CSIC-08-14 
Version 15th June 2010 

Separating expansion from contraction in spherically symmetric models with a 

perfect-fluid: 

Generalization of the Tolman-Oppenheimer-Volkoff condition and application to 

models with a cosmological constant 

Jose P. MimoscQ 

Departamento de Fisica, Faculdade de Ciencias da Universidade de Lisboa, Centra de Astronomia e Astrofisica, 
Universidade de Lisboa, Av. Gama Pinto 2, 1649-003 Lisboa, Portugal 

Morgan Le DeUiovJl 
Instituto de Fisica Teorica UAM/CSIC, 
Facultad de Ciencias, C-XI, 
Universidad Autdnoma de Madrid 
Cantoblanco, 28049 Madrid Spain 



a, Portugal 



Filipe C. Men£@ 

Centro de Matemdtica, Universidade do Minho, Campus de Gualtar, 4710-057 
(Dated: Received 2 November 2009; Accepted 5 May 2010) 

We investigate spherically symmetric perfect-fluid spacetimes and discuss the existence and sta- 
bility of a dividing shell separating expanding and collapsing regions. We perform a 3 + 1 split- 
ting and obtain gauge invariant conditions relating the intrinsic spatial curvature of the shells to 
the Misner-Sharp mass and to a function of the pressure that we introduce and that generalizes 
the Tolman-Oppenheimer-Volkoff equilibrium condition. We find that surfaces fulfilling those two 
conditions fit, locally, the requirements of a dividing shell and we argue that cosmological initial 
conditions should allow its global validity. We analyze the particular cases of the Lemaitre-Tolman- 
Bondi dust models with a cosmological constant as an example of a cold dark matter model with a 
cosmological constant (A-CDM) and its generalization to contain a central perfect-fluid core. These 
models provide simple, but physically interesting illustrations of our results. 

PACS numbers: 98.80.-k, 98.80.Cq, 98.80.Jk, 95.30.Sf , 04.40.Nr, 04.20.Jb 



I. INTRODUCTION 

Models of structure formation generally assume that 
small local iiihomogeneities grow due to the gravitational 
instability, so that the overdensities collapse and even- 
tually form the "bound" structures we observe in the 
present universe. Underlying this viewpoint is the idea 
that the collapse of the overdensities departs from the 
general expansion of the universe. This approach often 
relies on the idea that a small overdensity can be ap- 
proached as a closed patch in an otherwise spatially flat 
Friedmann universe, and it claims that BirkhofFs theo- 
rem justifies that, on the one hand, its evolution is in- 
dependent from the outside universe, and, on the other 
hand, that the behavior of the outside Friedman universe 
is immune to the collapse of the closed patch (see e.g. 
[l|-[3||)-The collapse of overdensities has been extensively 
studied and most works have been focused on the study 
of the formation both of small structure (astrophysical 
objects) and of large-scale structure as the outcome of 
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the growth of small perturbations in a cosmological con- 
text. The latter subject comprises the relativistic and 
Newtonian analysis of the evolution of the fluctuations 
(see e.g. 0-0) and the study of the subsequent ampli- 
fication of the growing modes into the nonlinear regime 
resorting to numerical methods (see e.g. |8l-lll|). In the 
present work we consider spherically symmetric, inhomo- 
geneous universes with pressure and study the question of 
whether there exists a dividing shell separating expand- 
ing and collapsing regions. Our goal bears a connection 
to the general problem of assessing the influence of global 
physics into the local physics (l2l.Tl3|. One aspect of this 
problem that has always attracted great interest is the 
endeavor to explain the local inertial phenomena in a 
Machian sense (see e.g. [H) and, in fact, Brans- 
Dicke theory |Ra - fl9| stems from this problem. 

Another related aspect has been the study of the influ- 
ence of cosmic expansion on local systems. Einstein and 
Straus [2Q] were the first to study this problem by con- 
structing a global solution that resulted from matching 
the spherically symmetric vacuum Schwarzschild solution 
to an expanding dust Friedmann-Lemaitre-Robertson- 
Walker (FLRW) exterior across a hypersurface preserv- 
ing the symmetry. Bonnor has made several investiga- 
tions along this line (see e.g. [21]). In particular, he 
copresented an exact solution representing a local dis- 
tribution of electrically counterpoised dust embedded in 
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an expanding universe with zero spatial curvature [22j, 
showing that the distribution participates in the expan- 
sion. Among the generalizations of this model are set- 
tings that keep the spherical symmetry but generalize the 
interior source fields by considering, for example, Vaidya 
(see [23j and references therein) or Lemaitre-Tolman- 
Bondi (LTB) spacetimes (see HH). On a different 
context, Herrera and co-workers Him have studied the 
"cracking" of compact objects in astrophysics using small 
anisotropic perturbations around spherically symmetric 
homogeneous fluids in equilibrium. The latter references 
are concerned with the existence of a shell where there is 
a change in the direction of the radial force acting on the 
particles of the shells. Whenever this happens one has 
a cracking situation, a concept introduced by Herrera 
in Ref. [29j. The approach of these works is somewhat 
complementary to ours because it is not the full evolution 
that is depicted there, but rather the effect on particles 
as a result from a departure from equilibrium. 

In this work we use a different approach from all the 
works described above. On one hand, by making use of 
a single coordinate patch, we do not have to handle the 
matching problem. On the other hand, our approach is 
not perturbative. We adopt the formalism that has re- 
cently been developed in a remarkable series of papers 
by Lasky and Lun u sing ge neralized Painleve-Gullstrand 
(GPG) coordinates [3242JI- We perform a 3 + 1 split- 
ting and obtain gauge invariant conditions relating not 
only the intrinsic spatial curvature of the shells to the 



Misner-Sharp massj [m] but also a function of the pres- 
sure that we introduce and that generalizes the Tolman- 
Oppenheimer-Volkoff (TOV) equilibrium condition. 

In particular, we consider that the existence of a spher- 
ical shell separating an expanding outer region from an 
inner region collapsing to the center of symmetry, de- 
pends essentially on two conditions. The first condition 
establishes that there is no matter exchange across the 
shell. The second condition establishes that the general- 
ized TOV equation is satisfied on that shell, and hence 
that this shell is in some sort of equilibrium. The differ- 
ence with respect to the original problem where the TOV 
equation was introduced for the first time is twofold. Our 
result does not rely on the assumption of a static equi- 
librium of the spherical distribution of matter, and con- 
sequently does not assume that all the internal spherical 
perfect-fluid spherical shells are constrained to satisfy the 
TOV equation. In our case the generalized TOV equation 
is just satisfied at the dividing shell. Besides, the gener- 
alized TOV function depends on the spatial 3-curvature 
in a more general way than the original TOV equation. 
Furthermore, we shall characterize the dividing shell with 
kinematic quantities that provide a gauge invariant for- 
mulation of the problem. 



In order to illustrate our results we will analyze some 
particular cases. The simplest example is provided by the 
well-known Lemaitre-Tolman-Bondi dust models with a 
cosmological constant that can be seen as an example 
of a A-CDM model. A preliminary presentation of this 
work can be found in [36(. As a second case we consider 
generalizations of the previous model to contain a cen- 
tral perfect-fluid core. These models provide simple, but 
physically interesting illustrations of our results. 

An outline of the paper is as follows: Section II The 
GPG-formalism of Lasky and Lun: 3 + 1 splitting and 
gauge invariants kinematical quantities. Section III Ex- 
istence of a shell separating contraction from expansion: 
general conditions. Section IV Particular examples: Sec- 
tion IV A A-CDM model (LTB with a cosmological con- 
stant). Section IV B Perfect- fluid core in a A-CDM 
model. Section V Discussion of our results. 

We shall use units such that 8irG = 1 = c, and the fol- 
lowing index convention: Greek indices a, (3,... = 1,2,3 
while latin indices a, 6, ... = 0, 1, 2, 3. 



II. 3 + 1 SPLITTING AND GAUGE 
INVARIANTS KINEMATICAL QUANTITIES 

In this section we set the basic equations that we shall 
subsequently need. For comparison, we follow closely 
the formalism used by Lasky and Lun [33], while slightly 
generalizing their derivations for the explicit presence of 
a cosmological constant A. 



A. Metric and ADM splitting 

We adopt the GPG coordinates of Ref. [33| and per- 
form an Arnowitt, Deser and Misner (ADM [37]) 3+1 
splitting [2H in which the spherically symmetric line el- 
ement assumes a perfect-fluid timelike normalized flow 
n a := — a\7 a t = [— a, 0,0, 0] (n a n a = —1), defin- 
ing with its lapse N — a and its radial shift vector 
V M = (/3,0,0), an evolution of the spatially curved 

three-metric 3 g Mt , = diag (^j^ , r 2 , r 2 sin 2 with time 



(dQ 2 



ds 2 



de 2 



~a{t,r) 2 dt 2 + 



1 



1 + E(t,r) 



(P(t,r)dt + dr)* 



r 2 dn 2 . (2.1) 



The 3+1 approach uses the projection operators along 
and orthogonal to the flow 

h ab 



N a b 



-n n b , 



--g ah + n a n b . 



(2.2) 



1 also referred to as ADM mass when considering the mass of the 
whole spatial hypersurface. 



where h ab is the 3-metric on the surface £ normal to 
the flow. Those projectors are also used for covariant 
derivatives: Along the flow, the proper time derivative of 
any tensor X ab d is 



cd 



e yab 
n A cd\ei 



(2.3) 
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and in the orthogonal 3-surface, each component is pro- 
jected with h 



(2.4) 



Then the covariant derivative of the flow, from its pro- 
jections, is defined as 

n a -b = N b c n a . c + n a . h = -n b h a + \Qh ab + a ab 

+ 0J ab , (2.5) 

where the projection trace, the expansion of the flow, is 
= ra a a , the rate of shear a ab is its symmetric trace- free 
part and its skew-symmetric part is the vorticity uj ab . 

For perfect-fluids we have the Raychaudhuri propaga- 
tion equation 



9 



* a a = - o© 2 + ri a h a - a ab a ab + Lo ab u ab 



-{p + 3P)+k. (2.6) 



where n — 87r. 

The quantity <d ab := ^C n h ab , where C n is the Lie 
derivative along the vector field n a , is the so-called ex- 
trinsic curvature and is given bj0 



6 afc =diag 
with N = 



0, 



l + E. 



P 



P 



atr" sin 



(2.7) 



2 l + E 
Its trace is the expansion scalai0 

(f3r 2 )' 1 L n E 



e = - 



which leads to the shear scalar 



2 1 + J5' 



1 r (P\ , 1 C n E 



Gl + E' 



3 a \r 



(2.8) 



(2.9) 



The 3-Ricci curvature tensor, which arises from fully 
projecting the Riemann tensor in accordance with 
Eq. (221, is 



Then, the 3-Ricci trace and trace-free 3-Ricci tensor de- 
rive from the 3-metric as 



! i? = -2 



(Br)' 



and 



6 r 2 



r (EY 
6 V r 2 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



where P% is diag [—2, 1, 1]. 

The trace and trace-free Hessian of a are given by 



(r 2 VTTEa'y (2.15) 



and 



-D^D v a - ^- 3 g^D c D c a = e{t,r)P^ 
a 3a 



with e = — - 



Yi+~e [YTTe y 



3a 



(2.16) 
(2.17) 



and where = h^V v is the notation for 3-covariant 
derivative used in Ref. (39j and in Ref. [l2j|. 

The Bianchi identity T b a a — can be projected along 
n b , giving 



n»T b a :a = - C nP - (p + P) 6 = 0. 



(2.18) 



while projections orthogonal to n b give the Euler equa- 
tion 



n a 1 b;c 




=> P' =-(p + P) 



B. The Einstein field equations 



(2.19) 
(2.20) 



3 P M „ = diag 



JiTEY^-\ E ' r - E ^ 

^E'r - e] sin ' ' 



(2.10) 



2 Recall that for a scalar C n = n a d a = —&t — @-dr\ |33|) called 
it K a i, but we preferred the Ellis convention for the extrinsic 
curvature. The prime denotes partial radial derivatives while 
the dot will denote from here on partial time derivatives. 

3 Note that we obtain a sign for and a different from that of 
Ref. [H. 



It is well known that the ADM approach separates the 
ten Einstein field equations (EFE) into four constraints 
and six evolution equations. Spherical symmetry reduces 
them to 2+2 equations. 

The Hamiltonian constraint reads, in the presence of a 
cosmological constant, 



5 i?+ |9 2 -6a 2 =1671-/9 + 2A, 



(2.21) 



the momentum constraint, restricted to the radial direc- 
tion by symmetry, 



(r 3 a)' = -y6' 



(2.22) 
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and the evolution equations can be reduced to3 



2C n Q - i 3 i? - 8 2 - 9a 2 + -D a D a a 
2 a 



= 24vrF-3A, (2.23) 
-£„a - a9 + e - g =0. (2.24) 



Using Eqs. (|2~%|) and (j2Uf in Eq. (l2~2"2"]l , one can simplify 
the latter into 



l + E a 2 "' 



(2.25) 



Using the guidance that, from Eqs. (|2.11j) and (|2.14l) . 
3 R + 12q eliminates derivatives in E, we can further sim- 
plify the combination of Eqs. [ ([2~25) | + 6 (|2~2"4"jl ] xr 2 /3 
with expressions from Eqs. (|2~%)> . ([2^5]) . (|2~TT|) . (|2~T4")l . 
and (j2TS"| as 



2r(l + E) (Ina)' -8irPr 2 + Ar 2 + 2r£ n 



a / \ d ( 
(2.26) 



Substitution of Eq. ([2~2l)j) into Eq. (|2~2"Tj) x r 2 /4 yields, 
together with Eqs. (ED, ([13]), (j2~TT|) . (j2"%)) . and 
r/2x (|2.26p . a Poisson-like equation that, integrated over 
r, defines a Misner-Sharp mass function |m| 

M' = 4vrpr 2 =*> M = 4tt / pa; 2 da; = r 2 (1 + P) (lna)' 
Jo 



4irPr 3 + -Ar 3 + r 2 L n I - ) , (2.27) 



which with Euler's Eq. (|2.20|) rewritten, for P ^ —p, 
leads to the expression 



M 

— + AirPr =C n I - ) + ^Ar 



«y 3 p + p 



P'. (2.28) 



The evolution Eq. (|2.26p can be recast to recognize the 
definition of ([2~2"T)l : 



E + 2 — + -Ar 2 = ' ^ 



M 1 
V + 3- 



(2.29) 



With Euler's Eq. ([2~2"0|) , the momentum Eq. ([2~2"S) be- 
comes 



(2.30) 



4 Note the sign differences in front of the Lie derivatives terms 
compared with [33t : our results give a sign for H which is con- 
sistent with the Raychaudhuri equation restricted to the FLRW 



while taking Eq. (|2.29p 's Lie derivative and using (|2.30[) 
with £ n ± = -f d r i = §/r 2 , then ^xEq. (I2~2~8l reads 



2/3 



C n M =AttPt z - 
a 



(2.31) 



Taking the positive (contracting) root of Eq. (|2.29p , the 
evolution Eqs. ax (12.31(1 and ax (|2.30p for M and E can 
be written in terms of time derivatives where we render 
explicit the Lie derivative (see footnote [5]) : 



M = a (M' + AnPr 2 ) \j 2— + i Ar 2 + E, (2.32) 



/ , 1 + E , \ / M 1 . n 

E = a { E+2 JTp~ p )VT + 3 Ar +E - (2 ' 33) 

This system is then closed with a choice of an equation 
of state (EoS). 



Generalized LTB 



Getting the metric (|2.ip into the generalized LTB 
(GLTB) form, as in [33] , requires a coordinate transform 
so that (3dt + dr oc dR. Taking t(T) and r(T,R), we have 
then the condition 

/3d T t + d T r = 0, (2.34) 

which becomes 

f3 = ~r. (2.35) 
Consequently, the line element (|2.1I) can be rewritten as 

„2 _ „ ,1m D \2 ro .\2 ,^2 i (5flr) 2 , „2 7 2 



ds^ = -a(T,RY(d T tydT z 



1 + E(T, R) 



dR z +r 2 dn 2 



(2.36) 



where E(T,R) > — 1 and we can freely absorb the time 
function in the new time by choosing t = T . Using now 
and ' for dx and Or, respectively, Eq. (|2.29p now reads 



=a ' I 2 — + -Ar 2 + E 



r 3 

and Eq. (j2~32"|l rewrites, using Eq. (g^g) 



(2.37) 



M =/347rPr 2 = 4nPr 2 a\ 2— + -Ar 2 + E, (2.38) 
V r 3 

while Eq. (|2.33p x r' rewrites 



Er > =2P 1 -±^P' = 2l±|p'a A /2^ + ^Ar 2 + £ 
p + P p + P V r 3 



and Euler's Eq. (|2.20l) xr' is unchanged, 

cV _ P' 
a _ p + P' 



(2.39) 



(2.40) 
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D. Remarks on A 

In all that precedes, the cosmological constant was kept 
explicit. However, from the EFEs, one can include its ef- 
fects in the total density and pressure as that of a fluid 
with pa = —-Pa = — • We then obtain expressions identi- 
cal to Lasky and Lun [33[ . It is interesting to note that 
the Misner-Sharp mass, in the explicit A formulation, is 
only referring to the initial, "A-less" mixture, while en- 
compassing the gravitational effects of the presence of A. 
From Eq. (|2.27p we can define the mass M tot and pressure 
term 47rP to tr 3 for the sum of the total perfect-fluid mix- 
ture plus A by taking Eq. (|2.27p for a perfect-fluid and 
setting A = 0. We can also interpret the sum of the to- 
tal mass and pressure terms as the mass of an equivalent 
dust model M e( j. We can then integrate the mass of A 
fluid and introduce the "Misner-Sharp mass" pressure 
term for the A fluid: 



+ 47rP tot r 3 =r 2 (1 + E) (In a)' + r 2 C n (^\ = M cd 

(2.41 



4-7T ,, A 



4TrP A r 3 



\Kr 3 . 



(2.42) 
(2.43) 



Thus we can rewrite the Misner-Sharp sum of the mass 
and pressure term from its components from Eq. (|2.27p : 



Ur 3 , 



M A + 4nP A r 3 = - ^Ar 3 + ^r 3 = -^Ar 3 



(2.44) 
(2.45) 



M + AnPr 3 =M tot + 4^P tot r 3 

1 a 3 A 3 1 
-Ar 3 H r 3 = 1 

2 6 3 

so M tot — M + M A and P tot = P + P A . In Sec. Hill 
unless stated otherwise, we will use M, p, and P to re- 
fer to the total values of the corresponding quantities, 
while we will adopt the notation M p f, p p f, and P p f to 
refer to the perfect-fluid quantities. We also wish to re- 
mark that although the mass evolution Eq. (|2.3ip refers 
to the "A-less" mixture mass and pressure, this conserva- 
tion equation holds for each component of a mixture of 
noncoupled fluids. We thus have for independent fluids 



M = M u 

fluid i 



(2.46) 
(2.47) 



fluid i 



C n M t =47rP,r 2 - = ±4 7 r J P J r 2 W2— + E. (2.48) 



M 
r 



III. GEOMETRICAL AND PHYSICAL 
CONDITIONS FOR THE EXISTENCE OF A 
DIVIDING SHELL 

In our spherical symmetric approach, we are looking 
for shells dividing expansion at all time from regions of 
mixed behavior involving periods of collapse. 



This leads to an investigation of the conditions for the 
dynamical separation of sections of matter trapped inside 
a dividing surface (physical condition). We will see that 
this approach is distinct from a purely kinematic sepa- 
ration of contraction from expansion (geometrical condi- 
tion) and will express the physical condition using kine- 
matic quantities. 



A. Misner-Sharp mass conservation 

In the previous section we have seen how the Misner- 
Sharp mass is evolving with the flow under Eq. (|2.31l) . 
We can thus define a surface for which this mass is con- 
served with respect to the flow: 



Vt, £ n M(t,r+(t)) =0 
M 

^>yt, E = -2 — , or P* = or r* = 0, 



(3.1) 



While the second case, P — 0, defines a dustlike layer in 
the perfect-fluid mix, and the third case, r = 0, is trivial, 
we shall concentrate on the first case. E = — 2— . In this 
case, from Eq. (|2.30p we get 



C n E 



r 



E- 



.1 + E 



P' = 0, 



(3.2) 



P + P 

so the shell is characterized by fixed curvature and 
Misner-Sharp mass. This implies that if a prescribed ini- 
tial P and p distribution is given such that there exists 
a shell where 



E± 



(3.3) 



then this shell can locally separate inner and outer re- 
gions that can be expanding and contracting differently. 
We call the separating shell a "limit shell," and denote 
it with *. In GPG coordinates the above condition is 



equivalent to ^ 
it to compute 



0, or to 



0. We can then use 



2M 

r* = a 

E 



C n M C n E 



2M 



r* = - 



M 

ClM 



M 



E 
E 



= 0, 



and 



a 



0. 



(3.4) 
(3.5) 

(3.6) 



so the limit shell appears as a "turnaround'!^ shell, in 
terms of areal radius. 

However, these conditions are coordinate dependent 
and give limited insight as to how they would express for 
different observers. This calls for a definition using gauge 
invariant quantities. 



See discussion in [lj Sec. 19, p. 77. 
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B. Expansion and shear 

Newtonian structure formation in spherical symmetry 
provides a natural limiting shell that is a locus separating 
at a given time expansion from collapse: the turnaround 
radius (see e.g. [40]). The definition of that locus is given 
by the vanishing of the expansion with respect to the 
flow. Nevertheless, this is not necessarily the case result- 
ing from condition 13. II Let us first start from the previ- 
ous mass flow definition and examine the corresponding 
expansion. 

In GPG coordinates [33], defining the flow by the 
shift/lapse vector, we can compute the expansion (the 
trace of the symmetric part of the projected covariant 
derivative of the flow vector), using Eqs. (|2.25p and (12.811 : 



2<?i 
a r 



At r* (for ^ = 0), we have nonzero expansion given by 

'P\' 



(3.8) 



The shear can also be expressed here from Eqs. (12. 9[) 
and (j2~2^| as 



1 

fl= 3 



£1 

a r 



(3.9) 



and we can then relate shear and expansion as [using 
Eq.r 



e, p r 

i I — + a \ — = L„r, 



(3.10) 



so on the limit shell, 

6* + 3a. =0 <S> (£„r) t = 0. (3.11) 



So, gTOV is equal to the radial acceleration or, more gen- 
erally, to the Lie derivative of ft /a, and hence Eq. (|3.15[) is 
the version in the GPG formalism of the classical Euler's 
equation of continuum mechanics. We also see that this 
gTOV acceleration relates to the force envisaged in the 
works of Herrera and collaborators [29l - l3lj | multiplied by 
(l + E)/(p+p), i.e., by (l-2Af/r*)/0> + p) at r = r*. 
We can then obtain local conditions that yield the TOV 
equation on the limit shell when 



sTOV. = & 



Cir =0 



e 



C n [-+a) =0. (3.16) 



We can further express gTOV in a form that reminds 
us of the FLRW Raychaudhuri equation by using (p) = 



(3 . 7) M/(47rr 3 /3), 



gTOV =J^P' + fr «p) + 3P) , (3.17) 
and for FLRW it reduces to 

gTOV FL =— r(p + 3P) = -r. (3.18) 



2. Dynamics of the limit shell 

We have seen that we could define the limit shell by 
only setting E+ = — 2M*/r* (so (5* = 0), so that 8* = 
3a*. Now, using Eqs. ([2^9"]) . (1232]) . (E22D, and (jXTBl . 
we find 



=p(^) +«gTOV (3.19) 

p =0 (p 1 - p— + - ) + a 2 gTOV, (3.20) 
\ a a J 



Generalizing TOV 



so on the limit shell, we have 



The TOV equation, following [H, emerges from 
Eq. (|2.28p in the static case. 

We now generalize the TOV equation by defining a 
functional gTOV from Eq. I|2~25|l as 



gTOV = 



1 + E P > 
Ppf+Ppf Pf 



M Pf 1 A 
— 7T Ar 



Using Eqs. \2AA\ and (|2.45|) we also have 
' /; /'' • \-l'r - 

r 



gTOV 



M 

72 



(3.12) 



(3.13) 



_P + P 

The definitions ([3TTD]) . ^TM) . and (jXT^I combine to yield 

2" 



gTOV 



C n \ — + a 



e 

3 +a 



-Ctr. 



(3.14) 
(3.15) 



(3.21) 
(3.22) 



^) =agTOV, 

=> P* =a 2 gTOV,. 

Recall that, in the LTB frame, p = — f, so this tells us 

r LTB ,* = - a 2 gTOV„ (3.23) 

and thus when gTOV^ = that shell has no acceleration 
and is therefore really static, as expressed in the origi- 
nal TOV equation. For completeness, we can reexpress 
Eq. HMD with Eqs. (|2.3ll2.30l3.13|) in GPG coordinates: 



2M 



rcPG, 



CiM ClE 



M 



gTOV, 



E 

. gTOV 2 - 



(3.24) 
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8r 2 dr 



From Eqs. (|3.7p we derive upon integration 

(3.25) 

where (/3/a) ro is a function only of t, which arises as the 
"constant" of the integration performed with respect to 
r, and which sets the value of (3 / a at r = tq. Using Eq. 
(l2~2"2"|) integrated directly, or Eq. (|3~2"5]) and (j3"TU|) , the 
latter result translates into 



a + f)=(7' '"' 



i 



(3.26) 



which is its gauge invariant expression. 
From Eq. (I3.25P we obtain 



r 



/3\ ^ 



Q 
1 

or 2 



3" \ T , 



r >To 



6 r 2 dr + 9 



>nW, ( ^ 

1 a 



r J ro 



er 2 dr 



= gTOV (3.27) 



This is the general equation that corresponds indeed to 
Eq. (21) of Di Prisco et al. [311, and it confirms their 
claim of a nonlocality of the radial acceleration. From 
Eq. (|3.25p we realize that this nonlocality also charac- 
terizes the radial expansion, as one should expect, and 
we further remark that a similar nonlocality is already 
present in the energy condition defining r* Eqs. (13. 3p and 
(|3~TTj) and in our gTOV condition Eqs. (|5TT3|| and ([535]) . 
since both implicate M, which is an integral between 
and r*. 

As j3 / a = r(0/3 + a) vanishes at r — r*, from 
Eq. (I3.25P one deduces that 



9 r 2 dr , 



(3.28) 



so that the integral on the right-hand side vanishes if 
the initial parameter rg (/3/a) vanishes at some interior 
value tq < tv This implies that 8 must vanish at some 
intermediate value tq < r < r*, since it has to change 
signs within the interval of integration. 

Likewise, when gTOV = 0, i.e. C n (/3/a), vanishes at 
some r, we derive from Eq. (|3.27p that 



r 2 | - 

a 



8 r 2 dr + 8 



1 



ro Jr <> 



„ 2Q t'P 



Qr 2 dr 



r -Jr 

So, at r = r*, the latter Eq. (|3.29p reduces to 



- / r -Jro 



Sr 2 dr 



(3.29) 



(3.30) 



and we conclude that the integral on the right-hand side 
vanishes if the term r^d^P / a) ro vanishes at 7'o. This 
result shows that the vanishing of the time derivative 
of the expansion thus occurs at least at one intermediate 
value between ro and r. In the case when dt(/3/a) ra rQ = 
at the center, we recover the result of Di Prisco et al. 
[3~D | , establishing the vanishing of the radial aceleration, 
i.e. 8 = 0, at some < r < r*. 



3. Raychaudhuri expansion evolution 



From Eqs. (j2~2"Tj) and (|2~2"3"]l . with A included as a fluid 
component, we have in the GPG frame, 



-2£ n 8 - -8 2 - 12a 2 + -D k D k a =8tt (p + 3P) . 
3 a 



(3.31) 



and on the limit shell, that reads 



--8* - 28 2 + -D k D k a± =8vr (p + 3P) , (3.32) 

a a 

showing that this shell can still be dynamic. Using the 
Euler Eq. ([2~20]l . the Hessian (|2~T5)) gives 



l D f Dja= l±lp' 
a 7 p + P 



E' 



2(ar 2 )' 



l+E 
1 +£ p , 



P + P 



(3.33) 



Thus Eq. P3T|) reads 



2 2/3 



4tt (p + 3P) - 



28 



P' 



3(3 



2 P' 



1 + V 
P + P 

' + E -P'. (3.34) 



r p + P J p + P 



Here, we can recognize the first term of TOV. On the 
limit shell the above equation reads 

1 • P' 



a -8,-8-4.(p + 3P)- 2(p + p) 



1 + P 
p + P 



P' 



P' 



l+E 



p + P J p + P 

(3.35) 



and we recast the Raychaudhuri equation for the FLRW 

case 



-£ n - ^- =4vr (p + 3P) 
= -3P-3P 2 . 



(3.36) 
(3.37) 
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4- Remarks on null expansion limit shells 

We now explore the consequences of having, in addition 
to (|3.1ip , the condition = Ofor the limit shell. In this 
case, the shear must also vanish on the shell and 



=0, 



(3.38) 



which constrains the gradient of the generalized velocity 
field p/a. 

In addition, and most importantly, the Raychaudhuri 
Eq. (|3.34p shows that an initially expansion-free dividing 
shell is not likely to remain so, and will drift radially. If 
we impose the vanishing of C n Q in Eq. (|3.31|) . we derive 



-D k D k a* =4tt (p + 3P)+ 



(3.39) 



which then translates into a thermodynamic condition 
on the second-order derivative of P, which should induce 
a very specific and ad hoc local equation of state of the 
perfect-fluid, namely 



1 + E 
p + P J 



P'\ = - 4tt (p + 3P) + 



P' 



P 



2(p + P 



P* 



Pi. (3.40) 



We conclude that the case of a static, expansion-free, 
limit shell is very restrictive: for example, in the sim- 
plest case, discussed below, of an inhomogeneous A-CDM 
model, Eq. (|3.40|) induces a restrictive equation of state 
P = — p/3 on the shell, which is verified neither by the 
dust component nor by the A fluid, whereas the limit 
shell in this case derives from a staticity condition (see 
Sec. EES). 



IV. APPLICATIONS TO SIMPLE MODELS 

We now will illustrate the behavior according to the 
limit shell of simple models. First we will see how it 
appears in a A-CDM model, that is, a Lemaitre-Tolman- 
Bondi dust model with a cosmological constant. We will 
then look at more general models including perfect-fluids. 



A. Overdensity in a A-CDM model 

In what follows we consider a A-LTB model which, be- 
sides the bare LTB case, is exactly solvable, the simplest 
perfect-fluid model with a cosmological context departing 
from LTB and which satisfies the conditions for the exis- 
tence of an asymptotically r-static dividing shell. Indeed, 
as stated in [33|], choosing P = leads to the usual LTB 



solutions. Setting P = in Eq. (|2~35f implied M = 0, 
and it is somewhat remarkable that this mass is still con- 
served for each shell in spite of the presence of A. A 
gives a homogeneous pressure, which in Eq. (|2.40l) gives 
a 1 = so we can redefine adT = dT* into the line ele- 



ment (|2.36p . and finally in Eq. (|2.39[) . assuming no shell 
crossing r' ^ 0. We are therefore left with Eq. (|2.37p in 
the classic LTB form, with 



=2- 



M 1 



Ar 2 + E. 



(4.1) 



Adding a cosmological constant modifies the mass def- 
inition but not the dust equation of motion. However, 
we have an extra term that leads to a different dynam- 
ics. We can thus write the Raychaudhuri-like equation 
corresponding to time derivation of Eq. (|4.1|) : 



M A 



(4.2) 



and this shows there exists a radius without acceleration 
for strictly positive A, contrary to pure dust. However, 
the first integral (|4.ip suffices for analysis of what hap- 
pens to each shell (with fixed R). 



1. Kinematic analysis 

The Friedmann-like Eq. (|4.ip can be used to get the 
dynamics in a purely kinematical way. It can be ex- 
pressed with a polynomial 



A. 
"37 



3E 



6M\ A 



(4.3) 



which roots (given in Appendix [SJ should obey the ef- 
fective potential equation 



2M A , 

E = V(r) = -r 2 . 

r 3 

Since r 2 > 0, we have the condition 
E >V(r). 



(4.4) 



(4.5) 



The motion of a given shell over time thus follows E = 
const curves above the effective potential V. Roots, 
the points of changing direction, translate as geomet- 
ric intersections between those curves and V. The ef- 
fective potential admits one real negative root (0 en- 
ergy/curvature) at 



6M 



(4.6) 



6 M can be understood as the mass of the dust alone but interact- 
ing with A, see Sec. Ill Dl 
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' 2M A„2 



E, 





Figure 1. Kinematic analysis for a given shell of constant 
M and E. Depending on E relative to £/u m , the fate of the 
shell is either to remain bound (E < < i?u m ) or to escape and 
cosmologically expand (E> > Ey lm ). There exists a critical 
behavior where the shell will forever expand, but within a 
finite, bound radius (E — Eum, r < Him)- The maximum 
occurs at mm ~ USM/A. 



Figure 2. Phase space of a shell of fixed M and E. The 
scales are set by the value of ru m = \/3A-//A while the actual 
kinematic of the shell is given by E. 



admitting only one real root; the acceleration is always 
positive for 



and one double solution at its horizontal tangent (V = 0) 



''in 



3M 
"X' 



for which the value of E becomes 

E iim =-(3Mf 3 A 1 /\ 



(4.7) 



(4- 



It can easily be shown that any shell standing at r% 
with En m will automatically be a limit shell 



2M< 



rib 



tot.lii 



M ■ 



At- 3 
6 lira 



E 



lim 



E ] 



lim 



3M 
'E~ 



lim 



(4.9) 



and calculating its gTOV, using the definition of 
Eq. (|3~T3)l and recognizing Eq. l[3~2]) . 



r z 3 



(4.10) 



that such a shell will be r-static (gTOV lim = —r\i m = 0). 
The effective potential analysis is shown in Fig. Q] 

We can thus reconstruct the phase space of that shell 
in the (r, r) plane. Above the energy E\i m , there is only 
one root in the negative region; thus the flow is qual- 
itatively defined by its initial conditions. At E\i m , the 
double positive root gives a repulsive point, thus a sad- 
dle, while, below Ey lmi the pair of roots give closed and 
open orbits as shown in Fig. [5J 

The Raychaudhuri-like equation can also be expressed 
with a polynomial 



A 

= 3r* 



3M\ _ A 



(4.11) 



r > 



3M 
"A" : 



(4.12) 



thus at infinity (cosmological constant dominates, and M 
is monotonous in r). Therefore, at this root, there exists 
a limit radius beyond which there is no recollapse: 



nim(-R) 



3M(R) 



A 



(4.13) 



Note that this radius corresponds to the saddle point, 
which initial energy radial profile is fixed with ini- 
tial conditions for the mass distribution E\i m (R) = 
- {3M(R)) 2/3 A 1 / 3 . Therefore the last intersection be- 
tween the initial curvature profile, set by combining ve- 
locity and mass profiles, and this saddle point profile 
yields a global shell beyond which there is no recollapse, 
recovering separation of expansion from collapse. Ex- 
plicit exact solutions for this ALTB evolution model are 
shown in Appendix [B] It is nevertheless crucial to realize 
that the selection of the limit shell from initial curvature 
does not entail necessarily that it should start as r-static. 
Indeed the opposite should be true in general, as can be 
seen in Eqs. (|4.1[) using E\ iln , R\ im in (|4.4p . and Fig.Q] for 
any choice of the initial i2u m < n™, the radial velocity 



Aim =^lim — V"(i?ii m ) > 0, 



(4.14) 



so it appears that the r-static behavior of the shell should 
only emerge asymptotically as it approaches zero veloc- 
ity for infinite time. The selected limit shell therefore 
agrees with the conditions (|3.11l3.16p only at infinity in 
time, and is traced back to initial conditions owing to the 
A+dust conservation of M and E in time. More general 
fluids should not always allow for this conservation on the 
limit shell; however, once a shell verifies Eqs. (13. 1113. 16p . 
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Figure 3. NFW with background density profile 



its staticity guarantees that it should verify it at time in- 
finity. It is remarkable that the existence of the limit shell 
only matters at time infinity, suggesting that a weaker 
definition than (|3.1ip and (13.161) . should be a sufficient 
condition. 



2. Time dependent TOV 



The shape of Eq. (|4.10l) shows that, at the root of 
the Raychaudhuri-like polynomial, gTOV = and that 
it is positive inside and negative outside. The trapped 
region is thus characterized by gTOV > 0. We can also 
compute, using M = A-k (p) r 3 /3, 



gTOV' 



47T p 



(4.15) 



so TOV is a decreasing function of r [for r' > 0, a fair 
assumption as seen when r(t = 0) = R], except in regions 
where p > | ((p) + pa), that is, in density peaks. It is 
also a time dependent function through the evolution of 
r: 



gTOV 



/2M 



2M A 



(4.16) 



and thus for a given shell, it increases with time for ingo- 
ing dust shells and decreases for outgoing ones. The main 
point is that with dust, turnaround shells have r-static 
gTOV, and that balanced shells (between their mass pull 
and that of A) verify the TOV equation and are thus 
static. 



3. Expansion and shear 

From the definition (|3.9p of the shear, we see that in 
the GLTB model under consideration 



1 (r 1 


r 


3 \r' 


r 



(4.17) 



Figure 4. NFW with background Eu m and an example of E 

" 10 and n=e 9 . 



profile given by Eq. (|4.22|l . for -E m i n = — 1 + e 



where we now denote by a prime the derivative with re- 
spect to the GLTB radial coordinate R (i.e., dr — dR/r 1 ). 
Using Eqs. (|4. 1 [) and (|4.2p we, then, derive 



1 



6^B +2 M + A r 2 



2E 
r 



(4.18) 



It is then possible to verify that this quantity does not 
vanish in general when r — > r*. It does vanish if the 
expansion O also vanishes at the locus where j3 /a = 0, 
i.e., at r = r*, as we have commented in Sec. IIIIB4I 



4- Examples of initial density 

It is obvious then that initial conditions are crucial to 
determine the existence of a separating shell in the ALTB 
model since they set the profile of E and that of Eu m . A 
single crossing of the two curves ensures locally the exis- 
tence of such a shell, while its global effect remains if the 
initial conditions do not foster shell crossing. This is the 
case if there is only one crossing from bound to unbound 
E of -En m . More complicated cases will be examined in 
a future work. We now proceed with examples of initial 
density profiles and then deduct the conditions on the 
corresponding curvature profile for a limit shell to exist. 

a. NFW with background: The choice of a so-called 
Navarro, Frenk and White (NFW) density profile [4l| is 
motivated by their prevalence in large cosmological dark 
matter haloes (Le Delliou [42j, and references therein). 
If we initialize the halo with such a density profile, with 
concentration 1/Rq and inflection density po/4, placed 
on a constant background pb, we can compute the corre- 
sponding mass profile. The density profile, as illustrated 
in Fig. [31 is given by [4l| 



Po 



_R (l + JL 

Ro \ T Ro 



Pb- 



(4.19) 
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The corresponding mass then reads 



M = 47T 



r oPo 



In 1 



R 
'o 



R 



R + r 



R } 



Pb- 



(4.20) 



Now armed with the expression for the maximum en- 
ergy function, the double root solution above, we can 
obtain from Eq. (|4.8p the bound upper limit for the ini- 
tial energy/curvature profile that separates between ever- 
expanding and bound shells 



)P0 



In 1 



R 

ro 

R 3 

- pb T 



R 



R + r 



2/3 



(4.21) 



Figure 0] shows that profile corresponding to the NFW 
with background mass. We then propose an example for 
the E(R) profile, motivated by its cosmological Fried- 
mann asymptotic curvature and its simple radial evolu- 
tion from bound to unbound, as 



E(R) = - 4£ n 



1 - 



.22) 



where r\ > and — 1 < E m i n < 0, chosen so that E 
crosses En m near its constant density region. With the 
asymptotic constant density and Friedmann negative cur- 



vature (E 



tR 1 



-fcoo-R 2 ), these initial conditions 



model well a collapsing structure in an open background 
of curvature radius The resulting curves are shown 
in Fig. [H We have here an example where shells with 
E < Bii m are trapped inside the limit shell defined by 
the intersection of the two profiles. Moreover, that limit 
shell in the case of dust with A has been shown to be 
static. Thus, with this set of physically motivated initial 
conditions, the limit shell defined in this way delimits a 
constant region of collapsing mass, separated from ex- 
panding shells. 

b. Cosmological background with power law overden- 
sity: The most natural cosmological initial condition is 
a power law overdensity, with or without cusp, upon a 
uniform background with an initial Hubble flow (Le Del- 
liou [111). The uniform background and initial Hubble 
flow ensures the asymptotic solution starts FLRW. In this 
second example of initial conditions, we explored both 
density profiles but illustrate only the cuspless case as 
it is more observationally sounded (Le Delliou (42|, and 
references therein). The density profiles, as illustrated 



for the second case in Fig. [SJ are given by (e > 0, and in 
the first case e < 3 for a finite central mass) 



R 

P=P ° life 



, , R 



Pb, 



Pb- 



(4.23) 
(4.24) 



tog(p(R)) 



Figure 5. Power law density profile without cusp and with 
background 
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Figure 6. Power law density without cusp + background 
inlog(— Sum) — log(_R) and log(— E) — log(_R) scales 



Observations of the cosmic microwave background would 
imply the choice of initial time at recombination and am- 
plitudes of the order of po ~ 10~ 5 /?fc (see Le Delliou [42l |. 
and references therein). The corresponding mass then 
reads, for the cuspy profile, 



M cusp = inrlpo 



H£) 



3-e 



e = 3 

47T , 

< e < 3 I 3 



(4.25) 



and for the profile with constant density in the center 
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The resulting boundary profile for E again follows 
Eq. (|4.8p . using the obtained mass profiles. Taking an 
initial Hubble flow, R = HiR, the E(R) profile is then 
defined by Eq. (I4.1J) to be 



E{R) = I Hf - ^ ) R 2 - 



2M 



(4.27) 



The resulting comparison between E and Eu m for the 
noncuspy case is shown in Fig. H3 Once again, the in- 
tersection defines a static limit shell for which ni m = 
and gTOV = 0, all shells inside it are in the 



2A/ tot 



kinematically bound region of Fig. [TJ while those outside 
are in the free region. Initial conditions ensure they will 
expand in a quasi-FLRW manner. 

These examples illustrate that cosmologically moti- 
vated initial conditions lead to a clear separation between 
expanding and collapsing regions. Therefore for these 
systems, expansion ignores the effects of collapse, and 
conversely the details of the collapsing region can ignore 
the presence of a background expanding universe. 



B. perfect-fluid core in a A-CDM model 

Before examining the possibility of existence for a limit 
shell inside a perfect-fluid in a sequel paper, where we 
shall present an ansatz for a perfect-fluid inhomogeneous 
core in a Friedmann environment, let us turn to the con- 
figuration where a perfect-fluid ball is surrounded by (a) 
vacuum with a cosmological constant, and (b) dust and 
A. 



1. Pure A exterior 

In the same way as [33j did for a perfect-fluid sur- 
rounded by a A = vacuum, we can examine the inter- 
face between the perfect-fluid and the A vacuum. In the 
latter region, both the pressure radial derivative P' = 
and the sum p\ + Pa = for all time and place by def- 
inition of A. In the same way as [33| showed for such 
a configuration with A = vacuum, such a simple in- 
terface implies, through Eqs. (|2.40[> and (|2.30[t . that the 
energy and lapse functions, E and a, are undefined there. 



lug(-E{R)) 

log(-E lim (R)y 



E > Eu 



tq escape to oc 
tq collapse to C 



if r j > 
if h < 



log(R) 



\E(S) > I) 



Figure 7. rn m < rg 1 < rg 2 case for a dust layer with A. Full 
space A-CDM diagram for log(— Eu m )— log(ii) and log(— E) — 
log(i?) in dashed line. This region is characterized by E > 
-Eiim, so the dynamical analysis of Fig. [T] yields continuation 
of initial velocities directions. 



These equations show that only if the fluid's pressure ra- 
dial derivative P' vanishes faster than p + P can E and 
a remain defined. This condition sets an unusual bound- 
ary constraint to the perfect-fluid's EoS (simple linear 
EoS do not agree with it), but it is more fruitful to point 
out that such behavior mimics that of a vanishingly thin 
layer of A-dust. Thus, the transition between the two 
regimes gives rise to an inescapable A-dust atmosphere, 
however vanishingly thin, as was found in the pure vac- 
uum case [33[ . We have two free boundaries, rg 1 (t) where 
the pressure vanishes and rg 2 (t) > rg 1 (t) where the den- 
sity vanishes, at which the EoS is defined as 



f(p,P) forre[0;r ai ] 
P for r e [r^jrgj . 



(4.28) 



Evolution of rg 1 (t) and rg 2 (t) follows from setting, re- 
spectively, P = 0, then P = p = in Eqs. ([2~3"2"]l . (E32J), 
and (|2.40p , to evolve those radii from initial conditions. 
The continuity of the curvature through both boundaries 
imposes again 



lim — lim 



{E(t,r)}=0, 



(4.29) 
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E(R) > 

Figure 8. rg 1 < riim < rg 2 case for a dust layer with A, A- 
CDM for log(-_Bi im )-log(_R) and Iog(--E) -log(.R) in dashed 
line. The region with E < i?ii m is trapped by its set of ef- 
fective potentials and will recollapse that with E > En m , so 
the dynamical analysis of Fig. [T] yields continuation of ini- 
tial velocities. The separating shell remains in between those 
regions. 



which can be used to transmit the value of the mass pa- 
rameter from the outer Schwarzschild-de Sitter spacetime 
down to the perfect-fluid boundary curvature. 



2. Limit shell 

At this stage, the possibility opens for a limit shell in 
the A-CDM atmosphere of the core, provided that such 
a shell verifies in conjunction Eqs. p.3[) . or equivalently 
(|3.1ip . and (|3.16p . which is only possible in a positively 
curved region. Given the surrounding Schwarzschild- 
de Sitter environment, the positive curvature require- 
ment is at least locally filled near the outer bound- 
ary. There the analysis of Sec. pV A[) applies fully to 
yield, given initial conditions, the location of the previ- 
ously discussed static virtual shell. Recall that in the 
Schwarzschild-de Sitter region, E — — 2Ml>2 — A r 2 while 

& ' r 3 

Sum = — (3Mg 2 ) 2 ^ 3 A 1 / 3 = est; however, the analysis 
only applies in the presence of dust, thus between rg 1 
and rg 2 . Owing to the preservation of continuity in M 
and E at rg 1 , whichever behavior the perfect-fluid may 
have, it will be confined by that of the previously ex- 
plored A-CDM at its boundary. 

Let us exhibit examples of such configurations: we 
can start from a similar example as presented in Sec. 
PV A 4[) . Nevertheless, to preserve curvature continu- 
ity (|4.29p . the initial velocity at rg 2 should go to 0, and 
therefore the previous E profile should be modified ac- 
cordingly. Then we are faced with three possibilities 
due to the location of the dust layer boundaries com- 
pared with the limit shell in the full space dust model: 
Him < r dl < rg 2 , r dl < r lim < rg 2 , or r dl < rg 2 < r lim . 
Those cases are illustrated, respectively, in Figs. [71 [5] and 
[9] In the first case, the dust layer locates above the maxi- 
mum of their effective potential (|4.4[) so their initial veloc- 
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Figure 9. rg 1 < rg 2 < rn m case for a dust layer with A, A- 
CDM for log(-_B Um )-log(i?) and log(-_B) -log(i?) in dashed 
line. This region is characterized by E < Eu m , so the dynam- 
ical analysis of Fig. [T] yields eventual recollapse. 



ities gives the direction of their unhindered asymptotic 
behavior; i.e. an initially expanding dust layer should 
expand forever. If a separating shell exists, it should lie 
within the perfect-fluid region. The second case shows 
the existence of a separating shell, the perfect-fluid being 
bound by the eventual recollapse of the rg t shell, while 
some of the dust shell will expand through the vacuum 
region and eventually squeeze it to infinity. In the third 
case all the dust shells locate below the maximum of their 
effective potential (|4.4p so the whole mass will eventually 
recollapse, as if the separating shell was virtually located 
in the vacuum region. 

Now sending the rg 2 boundary to infinity, we can ex- 
pand the dust layer accordingly and so long as Sec. 
PV A[) 's analysis yields a limit shell within the dust re- 
gion, the perfect-fluid shall be contained by the collapsing 
inner boundary (i.e. the third case disappear and we are 
left with cases rn m < rg 1 and rg 1 < rn m as treated in 
Figs. [JJ and HI). 

In this section we have found that the presence of a 
cosmological constant does not modify the need for a dust 
layer around a perfect-fluid core surrounded by vacuum. 
We have also given examples of limit shell separation 
behaviors for appropriately set initial conditions in the 
dust layer with A. We have even hinted at that possibility 
inside the perfect-fluid from the dust behavior, although 
such study should be left for a sequel paper. 



V. SUMMARY AND DISCUSSION 

In the present work we have considered spherically 
symmetric, inhomogeneous universes in order to ascer- 
tain under which conditions a dividing shell separating 
expanding and collapsing regions exists. This endeavor is 
important in relation with the present understanding of 
structure formation as the outcome of gravitational col- 
lapse of overdense patches within an overall expanding 
universe. 
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We have addressed this problematic by resorting to 
an ADM 3+1 splitting, utilizing the so-called gener- 
alized Painleve-Gullstrand coordinates as developed in 
Refs. f32l . [33l | . This enables us to follow a nonperturbative 
approach and to avoid having to consider the matching of 
the two regions with the contrasting behaviors [44|. We 
have found local conditions characterizing the existence 
of a dividing shell. We have related these conditions 
to a gauge invariant definition of the properties of the 
dividing shell. These require the vanishing of a linear 
combination of the expansion scalar and of the shear on 
the shell, as well as that of its flow derivative. In GPG 
coordinates, it summarizes as a vanishing of both first- 
and second-order flow derivatives of the areal radius. 

In order to illustrate our findings we have considered 
some simple examples of cosmological interest that pro- 
vide realizations of our results. We have considered a A- 
CDM model whereby we consider an LTB universe with 
dust and a cosmological constant. Notice that the simul- 
taneous consideration of the latter two components yields 
a perfect-fluid model for the combined matter content. 
Moreover it can be seen as a simplified model of a dust 
universe within a cosmological setting coarsely provided 
by A, which would then mimic the energy content of the 
background cosmological model with a rate of expansion 
much smaller than that of the pure dust collapse. 

We have chosen initial conditions motivated by cosmo- 
logical considerations and have discussed the existence 
of a dividing shell for those cases. We have also gener- 
alized a result of Ref. (author?) [33] for the case where 
a cosmological constant is present, which states that a 
perfect-fluid core embedded in a universe filled with a 
cosmological constant necessarily exhibits a dust transi- 
tion between the perfect-fluid inner region and the outer 
vacuum region. This permits one to envisage this case as 
a generalization of the former A-CDM examples. 

Finally we should mention that a thorough discussion 
of global conditions represents a much harder problem, 
and remains an open problem since this involves the full 
characterization of a partial differential equations prob- 
lem with boundary conditions in an open domain. 



ACKNOWLEDGMENTS 

The authors wish to thank Jose Fernando Pascual- 
Sanchez for bringing to their attention the work 
of the authors of Ref. [33j and for helpful discus- 
sions. The work of M.Le D. is supported by CSIC 
(Spain) under Contract No. JAEDoc072, with par- 
tial support from CICYT Project No. FPA2006- 
05807, at the IFT, Universidad Autonoma de Madrid, 
Spain, and was also supported by FCT (Portu- 
gal) under Grant No. SFRH/BPD/ 16630/2004, at 
the CFTC, Lisbon University, Portugal. F.C.M. 
is supported by CMAT, Univ. Minho, FCT 
Project No. PTDC/MAT/108921/2008 and Grant No. 
SFRH/BSAB/967/2010. J.P.M. also wishes to thank the 



provider of Grants No. PTDC/FIS/102742/2008 and No. 
CERN/FP/109381/2009. 



Appendix A: Roots of Psj(r) 

1. Roots for the polynomial 

The roots (ro) of Eq. (|4.3p proceed from the polynomial 
Psj- We change variables such that r = u + v and use 
the extra degree of freedom to choose to rewrite P 3 j = 
such that 



3M 



uv 
2 



E 

A' 
E 
A 



/3_M 



(Al) 
(A2) 



Solutions for the latter second degree polynomial come 
naturally as 



3M± + (3M) 2 
A 



-3M± yJSfc + (3M)' 
1 



i(277fe/3) 



(A3) 
(A4) 



We are left with six solutions for u and v, which are 
symmetrical and related by Eq. (|A1|) so uv being real, 
choosing u 3 as the positive square root solution, the cor- 
responding v 3 becomes the negative one while u and v 
are complex conjugate, so 



(3M ) 2 -f 



(3MY 



A 2 



E 
A' 



(A5) 



and therefore the roots are 



r*=o,±i = V-3M- 



-^J-3M - + (3M)V^ fe / 3 A /A 1/3 . 



(A6) 



2. Real root(s) 



For the positive discriminant, A = Mr- + (3M) 2 , there 
is only one real root for k = 0. A negative or null dis- 
criminant yields again the real k — root and two other 
real roots for k — ±1, since then v = u. We are then left 
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with the single real root, noting 




<=\I-M \/^ + (3M) 2 , (A8) 



m =- 



A 1 / 3 • 



(A9) 



and, when ^ + (3M) 2 < 0, the two other real roots 



a± = 



\ 



3M + iJ ( ^ - (3M) 2 (lTiiV3), (A10) 



^3M-i ^^ - (3M) 2 (l±h/3), (All) 

(A12) 



a± + a± 
2AV3 ' 



3. Signs of the real roots: 

So as to order the roots, it is necessary to look at their 
sign. This is important as r should be positive, r < 
being unphysical. Recall that M, A > 0, and E > 
-1. When A > 0, i.e. when E > - (3A/) 2/3 A 1 / 3 = 
-E'limi we have only one real root and tq > =>■ do > 



-oj. We always have -ag = y 3M + ^/-^ + (3M) > 

0. Supposing ao > (and thus > 0) then — ajpa 3 ) = 
^- > <^ £ > 0. Therefore, with the hypothesis E > 0, 
the condition r > implies ao > — <^> a 3 , > — a^ 3 
-3M > 3M. Hence for £>0we have r < 0. The same 
as for > E > — (3M) 2 ^ 3 A 1 / 3 , requesting ro > implies 
a > -Oq while -a^ > > a Q . Therefore, > E > 
— (3M) 2 ^ 3 A 1 / 3 always entails ro < 0, and we conclude 
that ro is always negative when E > — (3M ) 2 ^ 3 A 1 / 3 . 
The case when (3M) A < 1 is more interesting as we 
have three real roots for -1 < E < - (3M) 2/3 A 1 / 3 . Let 
us use the solutions of Eq. (|A6|) in the form 



u k +u k _ 2Re(w fe ) 



Ai 



We know that 



ul =-3M + i1 



(3M) 2 



(A14) 



so Im(w|) > and Re(n|) < 0. We can then 



rewrite ui 



pe ifk,3 w ith p 2 



and (f k ,3 G 



[ f + 2/ctt ; 7T + 2/ctt ] The values of Wfc are de- 



duced as Ufc 



p 3 e IVfc with 



<Pk.3 . 

3 ' 



Vft G 



[ I + 2fer . 1 + 2|tt j Each Ufc admits the same 

modulus, so the phases, each separated by 27r/3, give 
us the ranges and the order in which each root lies. The 
results are the following: 



^ g[| ; f]c[o ; f] 

<P+ G [ 7T - f ; 7T ] C [ § ; 7T ] 



99- G [ 







> r >0, (A15) 
r+ <0, (A16) 
r_ >0, (A17) 



and the order of the cosine (since r k involves the real 
part of Uk) yields — r + > ro > r_ > 0. This is agreeing 
with the analysis of Sec. IIV A II understanding that the 
negative root shifts from ro to r + through the A = 
point, and that below the horizontal tangent, ro is the 
exterior turning point while r_ gives the interior envelope 
of the effective potential. 

The above solutions give us then the explicit equations 
for the intersection of the effective potential with the cur- 
rent curvature involved in Eq. 14.11 



Appendix B: Exact solutions for an inhomogeneous 
ACDM 



The equation of motion admits analytical solutions in 
terms of hyperelliptic integrals (see also Lemaitre |43l|). 
From Eq. [gTTJ 



t-tn 



r U Er + 2M + f r 3 



dr; (Bl) 



however, in conformal time (dt = rdrf) 



A 



r' 2 =Er 2 +2Mr+ -r 4 , 



V - r)B 



R x jEr 2 + 2Mr+ fr 4 



-.dr 



zdr 



(B2) 



(B3) 



Given that the incomplete elliptic integral of the first 
kind is defined by 



(A13) F(x,k) 



dt 



y/(l - t 2 ) (1 - kH 2 ) JO 



X dt 



(B4) 



it is possible by a rational change of variable, z = ff^ 
to go from Pp to P4: 

P F (z(x)) = ((c - a)x + (d- b)) ((c + a)x + (d + b)) 
x ((c — ka)x + (d — kb)) 

x ((c + ka)x + (d + kb)) / (cx + d) 4 



{cx + d) 



(B5) 
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The solutions are therefore following, using cr + d 

f^H and dr = , ad - bc VJ dz: 

(a — cz) (a—cz) 



V~VB 



1 



1 



dr 



r JPAzj {cr + df 
(ad — be) 



(B6) 



We then just need to find a, b, c, d, k in terms of E, M, A. 
We already have the roots of P4 = P 3 jr |- from Appendix 
|A1 and we can write from Eq. (|B5[) 



d-b 



n 



r3 = 



c — a 
d-kb 

c~ ka 



-, r 2 



d + b 

c + a ' 
d + kb 

c + ka 



(B7) 



We can obtain expressions for d and 6, isolating them in 
the first and second pairs of roots: 



d= - 



ri(c — a) + r 2 (c + a) 



b = 



ri(c — a) — r 2 (c + a) 



r%(c — ka) + r±(c + ka) 
2 

7-3(0 — ka) — r^(c + ka) 
2k 



(B8) 



(B9) 



Equating the two ways of writing b + d, we obtain a linear 
relation between c and a, 

_ r 3 fc(l -k) + r 4 fc(l + k) - 2kr 2 
C - r 3{1 - k ) + 2 kr 2 ~r 4 (l + k) °' ^ 

Now recall that the factors of x 4 and x° in P4 are, re- 
spectively, 



( C 2 -a 2 )(c 2 -fcV) 
rir 2 r 3 r4 = 0. 



A 

3' 



(BH) 
(B12) 



The cosmological constant means from Eq. (|B11[) that 
neither c = ±a nor c = ±fca, while Eq. (IB12|1 entails 
that one of the roots is 0. If we choose r 4 — 0, then we 
have d = —kb and therefore, from Eqs. (|B8I) . d + kb = 
yields 



c _ n{l - k) - r 2 (l + k) 
a ~ n(l-fc)+r 2 (l + /c)' 



(B13) 



so with Eq. (|B10[) and r 4 = 0, we obtain a third degree 
polynomial in k (recall k ^ 1 for nondegeneracy of Pp) 

(fc - 1 )(f fc+ 2 ^-^3-r 2 r 3 y + 1 
V r i r 3 - r 2r 3 



2r 1 r 2 -r 1 r 3 -r 2 r 3 \ 2 l =Q 
rir 3 - r 2 r 3 / 



2rir 2 - nrg - r 2 r 3 / / 2rir 2 - nr 3 - r 2 r 3 ( 



r 2 r 3 - rir 3 



rir 3 - r 2 r 3 



(B15) 



We also can rewrite the condition (|B10[) to obtain a with 
Eq. (JBTTJ: the positivity of A in Eq. (|B"TT|l . 



A 



4k 2 (1 -fc 2 ) 2 (1 — A:) 2 r 3 + 4r 2 fc 



[r 3 - r 2 \ r 2 r 3 



[2r 2 k + (l-k) r 3 } 4 



(B16) 



imposes to choose r 3 > r 2 > 0, and thus 
a = ± [2r 2 k + (l-k) r 3 ] 



\ 



3[(1 — fc) 2 r 3 +4r 2 fc] [r 3 — r2]r2r 3 



2fc|l - k 2 \ 



(B17) 



We deduce then c from Eq. (|B10[) 
c = ±k[{l-k)r 3 -2r 2 ] 



\ 



3[(1 — fe) 2 r 3 +4r2*:] [r 3 — r 2 ]r 2 r 3 



2k\l - k 2 \ 



(B18) 



derive b from including the solutions (|B17IB10I) in its 
expression in Eq. (|B8l) 



b = T 



4r 2 k + (1 - k f r 3 n + [(l - k 2 ) r 3 ] r 2 



\ 



3[(1 — fe) 2 r 3 +4r 2 fe] [r 3 ~r 2 ]r 2 r 3 



2fc|l-fc 2 l 



(B19) 



and obtain d with our choice of r 4 = that induces 
d= -kb 



d = ± 



4r 2 k + (1 - k f r 3 n + fe [(1 - fc 2 ) r 3 ] r 2 



\ 



3[(1 — fe) 2 r 3 +4r 2 fc] [r 3 — r 2 ]r 2 r 3 



2fc|l-fc 2 | 



(B20) 



Inputting the values of the roots from Appendix [XI and 
the values of the transformation coefficients a, b, c, and 
d into Eq. (|B21[) yields the conformal time evolution so- 
lution that can be related to the cosmic time according 
to 



t-t B 



rdr)= r r d_ ( Fj^ky 

l/B Jr dr v ( ad - bc ) 



dr. (B21) 



Therefore there is an analytic solution to the ALTB 
model (see also Lemaitre [43]). 
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